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1.  Introduction 


The  motivation  for  this  study  is  to  arbitrarily  deliver  lethal  modules  to  prescribed  locations  on 
the  battlefield.  The  basic  idea  is  that  the  most  efficient  usage  of  energetic  material  may  be 
achieved  through  precise  arrangement  of  smaller  payloads  at  the  target  site.  The  launch  and 
lethal  mechanism  are  not  explicitly  considered  in  the  present  work. 

The  means  of  flying  these  small,  lethal  payloads  from  the  launcher  to  the  target  is  addressed  in 
this  study.  Of  particular  interest  are  concepts  that  reliably  produce  desired  terminal  patterns  at 
low  cost.  Feedback  control,  whereby  commands  to  control  the  system  are  based  on  feedback 
measurements,  is  a  classic  way  to  ensure  high  performance  in  a  complex  environment  such  as  a 
battlefield.  Indeed,  military  applications  have  driven  guidance,  navigation,  and  control  (GNC) 
technologies  for  aircraft1  and  missiles.2  5  Flocking  or  swarming  behaviors  of  unmanned  aircraft 
systems, 6-7  parafoils,8  and  spacecraft9  have  been  investigated.  Bennet  et  al.6  developed  a  safety- 
critical  guidance  algorithm  based  on  potential  fields  to  fly  aerial  robot  swarms  without  actuator 
saturation.  Collision-avoidance  and  formation-flight  algorithms  for  airdrop  of  supplies  from 
multiple  units  deployed  in  the  same  airspace  have  been  derived  and  verified  in  simulation  and 
experiment.8 

Guided  projectile  technologies  have  demonstrated10  12  substantial  progress  recently.  The 
emergence  and  subsequent  gun  hardening  of  Global  Positioning  System  (GPS)  technology  has 
enabled  precision  artillery  weapons.10  Fresconi11  used  GPS  feedback  in  a  nonlinear  model 
predictive  control  strategy  to  reduce  the  actuator  and  sensor  burden  on  a  rolling  artillery  and 
mortar  projectile.  Enhanced  maneuverability  while  maintaining  low-cost  actuation  components 
has  been  established  more  recently.12  These  efforts  have  largely  focused  on  a  single  flight  for 
larger  caliber  and  indirect-fire  applications. 

The  novel  contributions  of  the  present  work  include  formulation  of  a  unique  modular-delivery 
concept  and  derivation  of  guidance  and  flight-control  algorithms  for  this  concept.  A  parent-child 
delivery  strategy  is  introduced  to  meet  performance  goals  at  low  cost.  A  flight-dynamics-based 
control  algorithm  is  derived  and  optimization  theory  is  applied  for  the  guidance  law.  Models  of 
the  control  mechanism,  flight,  and  feedback  measurements  are  outlined.  These  models  are 
implemented  in  simulation  with  the  appropriate  input  data  along  with  the  guidance  and  flight- 
control  algorithms.  Results  illustrate  the  feasibility  and  key  features  of  the  novel  modular- 
delivery  concept  and  guidance  and  flight-control  techniques. 
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2.  Concept 


Options  for  the  launch  of  the  modular -delivery  concept  include  burst-fire  or  a  drop  from  an 
unmanned  aerial  system.  Regardless  of  the  specific  launch  mechanism,  this  work  assumes  that 
multiple  bodies  are  flying  simultaneously  in  a  relatively  close  proximity.  The  idea  of  a  parent- 
child  relationship  is  used  to  achieve  performance  metrics  with  simplified  components.  In  this 
paradigm,  a  parent  flight  body  features  higher  performance  components  (e.g.,  more  feedback 
sensors,  higher  precision  actuators,  faster  processors)  which  are  exercised  to  accurately  fly  to  the 
complex  target  and  gather  information  useful  for  the  child  flight  bodies.  The  children  are  shed 
from  the  parent  body  as  appropriate,  depending  on  the  mission,  and  are  equipped  with  simpler 
components  (e.g.,  reduced  sensor  suite,  minimal  actuation).  This  notion  exercises  swarming 
behaviors  to  enable  efficient  delivery  of  simple,  modular  bodies  to  desired  arbitrary  locations 
against  complex  targets. 

Consider  Fig.  1  as  an  example.  In  this  illustration  a  parent  body  contains  image-based 
navigational  technology  and  electromechanically  actuated  canards.  Image-based  navigation 
provides  extremely  accurate  targeting  for  the  parent  body  and  also  low  latency  information 
regarding  optimal  patterns  for  child  bodies  against  complex  target  arrays.  The  electromechanical- 
actuation  technology  allows  for  the  high  maneuverability  necessary  for  gliding,  loitering,  or 
intercepting  moving  and  defilade  targets.  The  benefits  of  these  components  are  often  offset  by 
added  cost  and  complexity.  As  a  means  of  producing  arbitrary  terminal  patterns  without 
exorbitant  cost,  the  child  bodies  in  Fig.  1  feature  a  simple  ranging  device,  such  as  a  radar  unit, 
and  a  ring  of  jet  thrusters.  As  the  body  rolls,  the  ranging  device  permits  the  relative  spacing  to 
surrounding  bodies  to  be  assessed.  Thruster  commands  can  then  be  issued  depending  on  the 
desired  pattern  of  parent  and  child  bodies.  Communication  between  bodies  (e.g.,  sensor  fusion, 
tasking  of  children  by  parent)  may  also  enhance  performance.  Trajectories  flown  to  a  complex 
target  array  using  this  parent-child  concept  are  shown  in  the  plot  on  the  right  of  Fig.  1. 
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Fig.  1  Illustration  of  parent-child  modular-delivery  concept 


3.  Flight  Models 


Flight  behaviors  of  parent  and  child  bodies  were  calculated  using  a  6-degree-of-fr eedom  model. 
Theoretical  details  of  the  body-fixed  model  are  provided  in  Murphy13  and  Fresconi  et  al.14 
Translational  and  rotational  kinematics  given  in  Eq.  1  and  2  relate  quantities  in  body-fixed 
coordinates  ([u  v  w]  and  [P  R  r])  to  inertial  coordinates  ([i  y  z]  and  [0  6  x p]) 

through  the  use  of  the  Euler  angles  ([0  9  0])): 
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The  translational  and  rotational  dynamics  presented  in  Eq.  3  and  4  determine  the  accelerations 
([ii  v  vv]  and  [p  q  f])  in  the  body-fixed  coordinate  system.  The  mass  properties  (mass  m, 
moment  of  inertia  tensor  /)  of  the  flight  body  are  required  along  with  the  aerodynamic  forces  and 
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moments.  The  aerodynamic  model  relating  states  of  the  airframe  (e.g.,  Mach  number,  angle  of 
attack)  to  aerodynamic  coefficient  data  (based  on  computational  fluid-dynamics  predictions  from 
Silton  et  al.15)  and  ultimately  forces  and  moments  is  not  provided  here.14  Gravity  is  also  used  to 
determine  the  overall  forces  and  moments  acting  on  each  flight  body  ([X  Y  Z ]  and 
[L  M  N ]): 
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The  control  forces  and  moments  for  the  child  body  are  supplied  by  a  continuously  variable 
lateral-thruster  array.  These  devices  use  some  means  (not  addressed  in  this  study)  to  smoothly 
throttle  the  thrust  produced  normal  to  the  spin  axis  of  the  flight  body  (i.e.,  no  control  of  axial 
component).  This  study  assumes  that  thrust  is  always  available.  Aerodynamic  canard  control 
with  a  first-order  model  for  actuator  dynamics  is  used  for  the  parent  body.14 

The  feedback  measurements  from  a  ranging  device  on  the  rolling  child  body  provide  the  location 
of  all  other  bodies  in  the  body-fixed  coordinate  system  of  that  specific  child.  The  3 -dimensional 
position  of  the  jth  body  with  respect  to  the  ith  body  in  the  ith  body-fixed  coordinates  can  be 
determined  in  the  6-degree-of-ffeedom  model  by  manipulating  the  inertial  positions  of  both 
bodies  and  the  Euler  angles  of  the  ith  body  according  to  Eq.  5.  The  parent  uses  a  seeker  and 
inertial-measurement  unit14  for  feedback  purposes: 
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The  1962  International  Standard  Atmosphere  model  provided  atmospheric  density  and  sound 
speed.  Winds  were  modeled  with  a  constant  magnitude  and  direction.  These  models  were 
implemented  in  simulation  and  integrated  using  a  fourth-order  Runge-Kutta  method.  Monte 
Carlo  capability  was  included  in  the  simulation  such  that  variation  or  uncertainty  in 
aerodynamics,  mass  properties,  atmospheric,  actuators,  feedback  measurements,  and  launch 
conditions  were  modeled.  Finally,  the  simulation  environment  was  made  extensible  to  an 
arbitrary  number  of  projectiles  in  flight  to  account  for  the  instantaneous  coupled  effects  of 
swarming  flight. 
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4.  Guidance  and  Flight-Control  Laws 


Algorithms  for  guiding  the  parent  were  presented  in  Fresconi  et  al.14  Here,  the  focus  is  on 
guidance  of  the  child.  The  relationship  between  feedback -measurement  output  and  control- 
mechanism  input  is  critical  to  flight-control  law  formulation.  For  this  reason  we  need  to  relate 
the  relative  position  from  the  ranging  device  to  the  thruster  force.  Fortunately,  Guidos  and 
Cooper16  have  provided  such  an  expression.  This  study  considered  the  flight  response  to  lateral 
impulses  similar  to  that  proposed  in  the  present  work.  The  equations  of  motion  governing  flight 
(Eq.  1-4)  were  simplified  for  this  situation  and  a  closed-form  solution  was  obtained.  A  version 
of  this  solution  is  shown  in  Eq.  6: 


Cnu  D  (xtgt  —  Et)  dF 
Cma  mV2  D  D 


(6) 


In  this  expression,  r  is  the  lateral  displacement,  CNa  is  the  aerodynamic  normal-force  slope 
coefficient,  C1Ua  is  the  aerodynamic  pitching-moment  slope  coefficient,  D  is  the  body  diameter, 
m  is  the  body  mass,  V  is  the  average  velocity,  xTGT  is  the  downrange  distance  to  the  target,  t  is 
the  time  of  flight,  dF  is  the  distance  of  the  thruster  force  from  the  body  center  of  gravity,  and  F  is 
the  thruster  force.  Thus,  Eq.  6  defines  the  lateral  displacement  (r)  of  a  body  flying  toward  a 
target  with  certain  characteristics  (CN  ,  Cm<x,  D,  m,  V,  xTGT,  t )  due  to  a  given  impulsive-control 

moment  ( dFF ). 


This  modular-delivery  concept  seeks  prescribed  distributions  of  child  bodies  around  the  parent. 
Using  Eq.  6,  a  predictive  control  strategy  can  be  built  around  control-force  inputs  based  on  the 
error  between  measured  and  desired  lateral  displacements.  Mass  properties,  aerodynamics, 
thruster  locations,  and  target  distance  can  be  obtained  prior  to  launch,  and  average  velocity  and 
time  of  flight  can  be  estimated  or  measured  in  flight.  These  values  can  be  collected  in  a  term,  AF, 
and  the  2  components  of  lateral-control  force  can  be  related  to  the  relative  position  feedback  as 
shown  in  Eq.  7 : 


y?j ]  =  _\Af  oip£j' 
z J  i  0  Af\  Fzb< 

.  IJ  J  L  J 


(7) 


Equation  7  provides  a  means  to  command  forces  on  a  given  child  using  ranging  data  to  one 
parent  or  child.  The  swarming  behavior  is  more  complex  because  a  method  must  be  devised  to 
optimize  force  commands  to  achieve  a  desired  pattern  from  multiple  bodies  surrounding  a  given 
child.  For  example,  a  desired  geometry  may  be  satisfied  between  the  child  and  the  parent  but  not 
between  that  child  and  another  child.  This  situation  is  depicted  in  Fig.  2  for  a  desired  circular 
pattern  of  5  children  (C1-C5)  equally  spaced  around  the  parent  (P). 


5 


Fig.  2  Resolution  of  desired  delivery  formation  for  multiple  bodies  in 
flight 

Optimization  theory  can  balance  the  desired  flight  formation  with  the  control  authority  available. 
These  problems  are  typically  cast  in  terms  of  a  cost  function,  J : 

J  =  eTQe  +  utRu  '  ^ 


Here,  the  control  error  is  e  and  the  control  action  is  it.  Weighting  matrices  for  the  control  error 
(Q)  and  control  action  (R)  permit  tuning  for  desired  flight-formation  performance  without 
prohibitively  high  control  commands. 

Next,  the  cost  function  is  minimized  with  respect  to  the  control  for  a  specific  formulation  and 
solved  to  obtain  an  expression  for  the  control: 
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For  the  present  problem  we  consider  only  the  control-error  portion  of  the  cost  function  shown  in 
Eq.  8  (i.e.,  R  =  0).  A  cost-function  variable  (e) ),  comprising  the  difference  between  the  lateral- 
displacement  feedback  (y;)  and  the  model  prediction  of  lateral  displacement  fromEq.  6  (Cu*), 
was  defined  to  solve  this  optimization  problem  for  the  ith  child  with  N  total  bodies: 
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Inserting  Eq.  10  into  Eq.  8  and  9  and  manipulating  yields  an  equation  for  the  control  of  the  ith 
child  as  given  in  Eq.  1 1 .  This  expression  essentially  follows  a  least-squares  problem: 
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The  final  step  is  to  resolve  the  lateral  commands  for  each  body  into  a  single  command  for  the 
control-force  magnitude  and  direction  of  the  ith  child.  This  could  be  done  through  schemes  such 
as  weighted  averaging;  however,  for  the  present  work  a  simple  average  was  used: 
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Figure  3  depicts  the  manner  in  which  control-force  magnitude  and  direction  are  resolved  into  the 
body-fixed  coordinate  system  for  each  child. 


Fig.  3  Control-force  magnitude  and 
direction  in  body-fixed 
coordinates  of  child 


5.  Results 


The  modular-delivery  concept,  guidance  and  flight-control  algorithms,  and  flight  models  were 
implemented  in  simulation  for  a  representative  40-mm  diameter,  fin-stabilized  airframe  with 
thrusters  located  3  calibers  ahead  of  the  center  of  gravity.  Swarms  were  launched  at  a  subsonic 
Mach  number  (150  m/s)  100  m  above  the  ground  against  a  target  400  m  downrange.  This 
situation  may  represent  a  launch  from  a  small  unmanned  aerial  system;  however,  the  intent  of 
this  study  is  to  demonstrate  the  utility  of  this  class  of  technology  across  a  wide  array  of 
applications. 

Numerical  experiments  were  conducted  to  assess  delivery  performance.  In  the  first  set  of  results, 
a  parent  and  single  child  were  flown  under  ideal  conditions  (e.g.,  no  errors  in  feedback 
measurements,  aerodynamics,  etc.)  with  a  desired  lateral  spacing  of  1  m.  Thus,  the  parent  should 
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impact  400  m  downrange  and  the  child  should  range  off  the  parent  body  yielding  an  impact 
400  m  downrange  and  1  m  off  to  the  side. 

The  flight  dynamics  of  the  parent  and  child  are  provided  in  Fig.  4.  Blue  curves  represent  the 
parent  and  purple  curves  are  for  the  child.  The  trajectory  as  a  function  of  downrange  distance  is 
contained  in  the  upper  plots  (cross  range  in  top-left  and  altitude  in  top-right).  The  cross-range 
flight  behavior  shows  that  the  parent  flies  along  the  line  of  fire  and  impacts  at  the  target.  Some 
oscillation  appears  in  the  child  cross-range  trajectory  and  the  impact  location  is  approximately 
1  m  to  the  side  of  the  parent,  as  desired.  The  trajectory  in  the  vertical  plane  is  almost  identical 
between  the  parent  and  child.  These  trajectory  results  suggest  that  this  guidance  and  flight- 
control  laws  are  successful  in  commanding  maneuvers  in  the  child  based  on  ranging 
measurements  off  the  parent  for  arbitrary  modular-delivery  patterns. 
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500 


200  300 
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Fig.  4  Flight  dynamics  of  parent  and  one  child 

The  pitch  (a,  solid  line)  and  yaw  (/ 3 ,  dashed  line)  angle-of-attack  histories  are  provided  in  the 
bottom-left  plot  of  Fig.  4.  Small  total  angle  of  attack  is  evident  in  the  parent  flight;  there  is  a 
positive  pitch  angle  of  attack  perhaps  because  the  target  is  slightly  beyond  the  ballistic  range. 
Interestingly,  the  child  pitch  angle  of  attack  follows  the  parent  history  with  small  oscillation, 
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which  is  consistent  with  the  vertical-trajectory  results.  The  yaw  angle  of  attack  for  the  child 
features  damped  oscillation  (as  seen  in  the  horizontal  trajectory)  to  some  appreciable  magnitudes 
(>5°).  The  Mach-number  histories  are  shown  in  the  bottom-right  of  Fig.  4.  The  higher  total  angle 
of  attack  of  the  child  yields  slightly  higher  velocity  decay. 

The  control  input  and  output  for  the  parent  with  one  child  are  provided  in  Fig.  5.  Lateral 
displacements  of  the  child  from  the  parent  obtained  from  the  ranging  device  on  the  child  are 
given  in  the  left-most  plot.  The  y-direction  (solid  line)  mimics  the  horizontal  trajectory  as 
expected  and  demonstrates  suitable  control  performance  (1-m  offset  at  impact).  The  z-direction 
(dashed  line)  data  indicate  small  displacements  as  desired.  These  feedback  measurements  are 
used  to  form  the  control-commands  input  to  the  system.  The  right-most  plot  has  dual  y  axes.  The 
left  y  axis  (blue)  is  for  the  maneuver  direction  ( <fic )  and  the  right  y  axis  (green)  is  for  the 
maneuver  magnitude  (F).  The  child  is  commanded  to  move  mainly  left  and  right  for  the  majority 
of  the  flight  and  demands  some  up  maneuvers  toward  the  end  of  the  flight  to  keep  from  falling 
short  of  the  parent  since  the  child  yaw-drag  component  was  higher.  The  magnitude  fluctuates 
when  the  child  flies  through  nodes  of  the  desired  set  point  but  overall  is  small  (<2N). 


Fig.  5  Feedback  output  and  control  input  of  parent  and  one  child 

Some  additional  control  metrics  are  detailed  in  Fig.  6.  The  upper-left  plot  is  the  time-varying 
control  gain,  AF.  These  data  show  how  a  higher  control  demand  is  called  for  earlier  in  the  flight 
regardless  of  the  error  between  the  feedback  and  model  prediction.  The  error  between  the 
feedback  and  model  prediction  in  the  y  and  z  directions  is  contained  in  the  upper-right  plot. 
Ideally,  these  errors  would  rapidly  converge  to  zero.  The  data  in  the  plot  illustrate  reasonable 
performance.  Oscillation  is  likely  due  to  using  a  position-feedback  measurement  with  a  force  (or 
acceleration)  control;  essentially  there  is  a  lag  between  when  measurements  are  collected  and 
when  the  response  appears  in  the  system.  Additionally,  the  guidance  law  is  framed  around  the 
idea  that  we  want  the  current  lateral  displacement  (from  the  feedback)  to  equal  the  final  desired 
lateral  displacement  (from  the  model  prediction).  These  oscillatory  effects  could  be  removed 
further  through  a  derivative  or  rate-feedback  term  in  the  controller  or  also  by  using  an  adjustable 
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finite  horizon;  for  example,  putting  the  “final”  downrange  distance  used  in  the  controller  (Eq.  6) 
closer  to  the  launcher  and  moving  this  farther  toward  the  true  target  at  each  control  update.  Even 
without  these  potential  enhancements  the  current  controller  performs  well  as  evidenced  by 
meeting  the  goals  of  the  desired  impact  patterns.  The  bottom  plot  gives  the  control  magnitude 
commands  in  the  y  and  z  directions  that  result  from  the  error  signals  in  the  upper-right  plot.  As 
expected,  larger  errors  in  the  y  direction  yield  larger  control  forces.  The  maneuver  magnitude  in 
the  right-most  plot  of  Fig.  5  is  the  root-sum-square  of  these  values  in  the  bottom  plot  of  Fig.  6. 


Fig.  6  Control  metrics  of  parent  and  one  child 

The  robustness  of  this  control  system  to  system  error  was  addressed  through  Monte  Carlo 
simulation.  Variation  or  uncertainty  in  aerodynamics,  mass  properties,  atmospheric,  actuators, 
feedback  measurements,  and  launch  conditions  were  modeled  as  described  previously.  A 
sampling  of  Monte  Carlo  results  for  a  parent  and  2  children  are  given  in  Fig.  7-9  using  a  format 
similar  to  that  in  Fig.  4-6  for  a  parent  and  one  child.  The  control  goal  for  this  situation  is  to  fly 
the  2  children  to  impact  1  m  to  the  left  and  right  of  the  parent. 

Figure  7  presents  the  flight-dynamic  behaviors.  Blue  curves  are  for  the  parent,  purple  curves  are 
for  the  first  child,  and  now  green  curves  are  for  the  second  child.  The  cross-range  data  in  the 
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upper-left  plot  are  interesting  in  that  the  parent  flies  off  the  line  of  fire,  perhaps  due  to  wind  or 
errors  in  the  feedback,  and  both  children  remain  offset  from  the  parent  by  1  m  to  either  side.  The 
control  metrics  are  also  met  in  the  vertical  trajectory  as  the  data  for  the  parent  and  child  bodies 
are  very  similar.  This  assessment  has  been  conducted  despite  the  difficulties  associated  with 
understanding  the  behavior  of  this  modular-delivery  concept  (e.g.,  the  need  to  simulate  the 
instantaneous  coupling  of  the  bodies  in  flight  with  high  fidelity).  These  trajectories  indicate  that 
the  modular-delivery  concept  and  guidance  and  flight-control  laws  are  robust  to  system  error. 


Fig.  7  Flight  dynamics  of  parent  and  2  children  with  uncertainty 


The  angle-of-attack  (bottom-left)  and  Mach  (bottom-right)  histories  in  Fig.  7  are  similar  to  those 
presented  in  Fig.  4.  The  yaw  angle  of  attack  for  the  2  children  is  mirrored  since  one  child  flies  on 
the  left  of  the  parent  and  the  other  is  to  the  right.  Mach  number  is  similar  for  both  children  since 
the  total  angle  of  attack  is  similar. 


The  control  input  and  output  are  given  in  Fig.  8  for  a  parent  and  2  children  with  system 
uncertainty.  Lateral-displacement  feedback  (left-most  plot)  is  similar  to  that  discussed  in  Fig.  5 
except  now  the  second-child  data  show  up  on  the  opposite  side  of  the  parent  from  the  first  child. 
The  control  commands  are  shown  in  the  right-most  plot  just  for  the  first  child.  These  data  closely 


11 


mimic  data  shown  for  the  parent  and  one  child.  The  only  difference  is  due  to  the  system  errors 
considered  in  the  present  Monte  Carlo  result. 


Fig.  8  Feedback  output  and  control  input  of  parent  and  2  children  with  uncertainty 

The  control  metrics  of  the  control  gain,  control  error,  and  control-force  magnitude  are  shown  in 
Fig.  9.  Again,  results  are  similar  to  those  for  the  parent  and  one  child.  Control  gain  linearly 
decreases  with  time  of  flight.  The  errors  and  control  forces  decrease  and  oscillate  with  time.  The 
majority  of  the  control  action  is  in  the  horizontal  plane  and  the  responses  of  the  2  children  mirror 
each  other. 

The  collective  flight  behavior  of  10  bodies  was  investigated  to  assess  the  extensibility  of  the 
concept,  guidance  and  flight -control  algorithms,  and  simulations.  One  parent  and  9  children  were 
flown  without  uncertainty  and  with  a  desired  pattern  of  1  -m  radius  from  the  parent  with  equal 
angular  (i.e.,  40°)  increments.  Format  for  presentation  of  these  results  follows  that  outlined 
previously. 
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Maneuver  Magnitude  (N) 


Fig.  9  Control  metrics  of  parent  and  2  children  with  uncertainty 

Flight  behaviors  are  collected  in  Fig.  10.  A  blue  line  denotes  the  parent  and  a  variety  of  colored 
lines  are  for  the  9  children.  The  trajectory  results  demonstrate  excellent  control  performance  as 
the  swarm  impacts  at  the  target  with  the  desired  pattern.  The  cross-range  and  angle-of-attack  data 
highlight  the  complexity  associated  with  controlling  the  flight  of  multiple  vehicles.  Each  child  is 
steered  through  pitch  and  yaw  angles  of  attack  to  yield  a  specific  lateral  displacement  at  impact 
(i.e.,  larger  yaw  angle  of  attack  for  larger  desired  cross-range  impact). 
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Fig.  10  Flight  dynamics  of  parent  and  9  children 

The  feedback  in  the  left-most  plot  of  Fig.  1 1  also  demonstrates  how  the  y  and  z  components 
scale  with  the  desired  pattern.  Some  oscillation  is  evident  in  these  data  and  the  feedback  for  all  9 
children  stabilizes  well  to  nearly  a  constant  value  for  the  latter  half  of  the  flights.  The  control 
commands  for  the  first  child,  shown  in  the  right-most  plot  of  Fig.  11,  replicate  the  data  for  a 
parent  and  one  child  in  Fig.  5. 
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Fig.  11  Left,  feedback  output  and  control  input  of  parent  and  9  children;  right,  control  commands  for  first 
child. 

The  control  metrics  are  provided  in  Fig.  12  for  completeness.  Inspection  of  these  data  shows  that 
the  behavior  of  these  data  for  10  bodies  follows  the  discussion  in  Fig.  6  and  9  for  fewer  bodies 
and  with  uncertainty. 


Fig.  12  Control  metrics  of  parent  and  9  children 
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6.  Conclusions 


Precision  delivery  of  smaller,  lethal  modules  on  the  battlefield  for  optimal  efficiency  could 
significantly  enhance  future  weapon  effectiveness.  A  major  challenge  to  realizing  this  vision  is 
the  means  of  reliably  achieving  prescribed  impact  patterns  without  excessive  cost.  A  parent- 
child  concept,  whereby  a  parent  body  possessing  higher-performance  components  guides 
multiple  child  bodies  equipped  with  simpler  components  to  the  target,  was  introduced. 

Novel  guidance  and  flight-control  laws  for  this  parent-child  concept  were  derived.  These 
algorithms  were  unique  in  using  the  flight  dynamics  to  appropriately  define  the  relationship 
between  control  input  (thruster-force  vector)  and  output  (lateral  displacement  from  other  bodies) 
and  optimization  theory  to  resolve  conflicts  in  formation  flight  for  the  child. 

High-fidelity  flight  models  were  utilized  to  successfully  demonstrate  the  concept  and  guidance 
and  flight-control  algorithms.  Capturing  nonlinearities  in  flight,  feedback  measurements,  and 
control  mechanisms  was  critical  to  proper  assessment  because  this  concept  instantaneously 
couples  the  motions  of  multiple  bodies  in  flight.  Simulations  (including  Monte  Carlo  analysis) 
were  challenging  from  a  computing  aspect  but  necessary  to  determine  robustness  to  system 
uncertainties  and  extensibility  to  larger  numbers  of  bodies  flying  to  arbitrarily  prescribed 
locations. 

Results  for  a  representative  system  demonstrated  the  feasibility  of  this  concept.  Despite  the 
challenges  of  accommodating  system  nonlinearities  with  reduced  feedback  and  simple  thruster 
control,  arbitrary  delivery  patterns  were  achieved  even  in  the  presence  of  uncertainties.  The 
children  flew  at  low  total  angle  of  attack  (no  evidence  of  flight  instabilities)  producing  a  slight 
decrease  in  Mach  number  from  the  parent.  The  lateral-displacement  feedback  from  the  ranging 
device  was  sufficient  for  control  purposes.  The  control  commands  featured  low  force  demands 
and  reasonable  directions.  Varying  control  gains  based  on  the  flight  dynamics  was  a  major  factor 
in  the  success  of  this  concept.  The  control  errors  converged  to  zero.  The  means  of  reducing 
oscillation  in  the  control  errors  and  response  include  adding  adjustable  finite  horizons  or 
derivative  terms  to  the  controller. 
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